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Abstract 

In quantizing gravity based on stochastic quantization method, the stochastic time 
plays a role of the proper time. We study 2D and 4D Euclidean quantum gravity in this 
context. By applying stochastic quantization method to real symmetric matrix models, 
it is shown that the stochastic process defined by the Langevin equation in loop space 
describes the time evolution of the non-orientable loops which defines non-orient able 2D 
surfaces. The corresponding Fokker-Planck hamiltonian deduces a non-orientable string 
field theory at the continuum limit. The strategy, which we have learned in the example 
of 2D quantum gravity, is applied to 4D case. Especially, the Langevin equation for the 
stochastic process of 3-geometries is proposed to describe the ( Euclidean ) time evolution 
in 4D quantum gravity with Ashtekar's canonical variables. We present it in both lattice 
regularized version and the naive continuum limit. 
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Introduction 

String field theory flU is believed to be the most promising approach to investigate non- 
perturbative effect in string theories. Recently, non-critical string field theories have been 
proposed for c = and for < c < 1 § §. Among these, some 

are based on the transfer-matrix formalism || in dynamical triangulation of random 
surfaces || . While some || are derived by using stochastic quantization method |10 



In the approach by stochastic quantization of matrix models, one can interpret the 
stochastic ( fictitious ) time as a Euclidean time coordinate in 2D quantum gravity. There 
is also another example of this particular observation, the proper time interpretation of 
the stochastic time, which was found in the course to study QCD^ in terms of the Nicoli- 
Langevin maps |]TTJ, that stochastic quantization of 3D Chern-Simons theory recovers 
the time evolution in 4D Euclidean Yang-Mills theory fli~2 |. These facts motivate us to 
interpret the stochastic time as the time coordinate in Euclidean 4D quantum gravity [|14 



In this paper, we illustrate how to apply stochastic quantization method to real sym- 



metric matrix models |T3[ and show that it leads to a field theory of non-orientable 
(non-critical) strings [EJ. The stochastic process defined by the Langevin equation in 
loop space describes the time evolution of the non-orientable loops on non-orientable 2D 
surfaces. The corresponding Fokker-Planck hamiltonian is a loop space hamiltonian of 
non-orientable string field theories. At the equilibrium limit, it deduces the Virasoro 
constraint equation for the probability distribution functional. The continuum limit of 
the field theory of discretized non-orientable loops is taken for the simplest one-matrix 
case ( c = ) and deduces the continuum field theory of non-orientable strings. 

Then we apply the strategy we have learned in 2D case to 4D Euclidean quantum grav- 



ity ||14|| . The Langevin equation for 3-geometries is proposed in the Ashtekar's canonical 
variables to describe the time evolution in 4D Euclidean quantum gravity in the sense, 
that the corresponding Fokker-Planck hamiltonian recovers the hamiltonian of 4D quan- 
tum gravity exactly. The stochastic time corresponds to Euclidean time in the temporal 
gauge, N — 1 and N l = 0. In this context, 4D quantum gravity is understood as a 
stochastic process where the scale of the fluctuation of "triad "is characterized by the 
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curvature at one unit time step before. The lattice regularization of the approach in 4D 
Euclidean spacetime is also presented to play the same game as what we have done in 
2D case with matrix models. 



Stochastic Quantization of Real Symmetric Matrix Models in Loop Space 

Let us start with the Langevin equation for one matrix model, 

AJWii(r) = -^£(M)<,(t)At + A&(t), 
S(M) = - -^N- a ' 2 tiM a+2 , 

(1) 

Mij denotes a real symmetric matrix. The stochastic time r is discretized with the unit 
time step At. We consider the discretized version of time evolution Mij(r + At) = 
Mij{r) + AMy(r) , with the Langevin equation for convenience of stochastic calculus 
and for understanding the corresponding stochastic process precisely. The discretized 
stochastic time development with At precisely corresponds to the one step deformation 
in dynamical triangulation in random surfaces. In the following argument, the specific 
form of the action of the matrix model is not relevant. The correlation of the white noise 
A£ij is defined by 

< A&(t)A6h(t) > 5 = Ar(8uS jk + fofy) . (2) 

It is uniquely determined^ from the requirement that (1) is transformed covariantly pre- 
serving the white noise correlation (2) invariant under the transformation M — > UMU^ 1 , 
where U denotes orthogonal matrices for the real symmetric matrix models. 

The basic field variables are loop variables <p n = tr(M n )A^ _1 ~t . Following to Ito's 
stochastic calculus WB, we evaluate 



n = <f> n (r + At) - n (r) , 

1 n-2 

= nXx{AMM n - x )N- 1 -^ + -n ^ u\ AMM h AMM" l ' '-j.V ; ^ + 0(Ar 3/2 ) . 

2 k=0 

(3) 

1 For an hermitian matrix Mjj in (1), the nose correlation is < A^j,-(r)A^/-;(r) >£= 2ATSnSjk ■ 



The terms in R.H.S. should be of the order Ar, thus we obtain 

n n ~ 2 1 
A0 n = Ar-{ <Pk<Pn-k-2 + (n - 1)— 0„- 2 } + Ar n ^ g a (pn+a + Mn-l , 

1 k=0 iV a=0 

ACn-i = ntr(A£M n-1 )7V _1- t . 

(4) 

The correlation of the new noise variables appeared in (4) is given by 
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< AC m _i(r)AC„_i(r) > € = Ar-nm < m+n _ 2 (r) > 5 , (5) 

The new noise is not a simple white noise but includes the value of the loop variable itself. 
In a practical sense, it might be tedious to generate the noise variable. We notice that 
<t>m+n-2(j) in R.H.S. of eq.(5) does not include the white noise at r but the series of noises 
up to the one step (stochastic time unit Ar) before. This means that the expectation 
value in R.H.S. should be defined with respect to the white noise correlation up to the 
stochastic time r — Ar. 

We also notice < A( n (r) >^= by means of Ito's stochastic calculus. In the context of 
SQM approach, the property of the noise yields the Schwinger-Dyson equation by assum- 
ing the existence of the equilibrium limit at the infinite stochastic time, or equivalently, 
lim^oo < A0 n (r) > ? = . We have, 

n n ~ 2 1 1 

<Pk<Pn-k-2 + 7;( n - !)T70n-2 + U ^ # a 0n+a >£= . (6) 
1 k=0 1 iV a=0 

The order of the noise correlation (5), l/N 2 , realizes the factorization condition in 
the large N limit. Therefore we obtain the S-D equation at large N limit for discretized 
non-orientable strings. 

j n-2 

7y < >£< 4>n-k-2 >5 9a < 4>n+a >£= . (7) 

1 k=0 o=0 

This shows that the S-D equation for non-orientable strings takes the same form as that 
for orientable strings at large N limit. The correspondence at the large N limit is exact 
if we define the corresponding hermitian matrix model by replacing all the couplings , 
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g a — > 2g a in (1). As a consequence, the disc amplitude in non-orientable strings is exactly 
the same as that in orientable strings. 

The geometrical meaning of the stochastic process described by the Langevin equation 
(4) is the following. The one step stochastic time evolution of a discretized loop, <f> n (r) — > 
0n( r ) + A0 n (r) , generates the splitting of the loop into two smaller pieces, 0& and n _fc_ 2 . 
The process is described by the first term in R.H.S. of (4). In a field theoretical sense, 
it is interpreted as the annihilation of the loop n and the simultaneous pair creation of 
loops, 0^ and <fi n -k-2- The first term in R.H.S. of (4) preserves the orientation of these 
loops, while the second term, which is the characteristic term of the order of for non- 
orientable interaction, does not preserve the orientation. Since the new noise variables 
in (5), A£„_i's, are translated to "annihilation" operators in the corresponding Fokker- 
Planck hamiltonian, the factor 2 in the correlation (5) for the new noise variables comes 
from the sum of the orientation preserving and non-preserving merging interactions. 
Namely, (5) describes the simultaneous annihilation of two loops m and 0„ and the 
creation of a loop 4> m +n-2- The geometrical picture allows us to identify the power "ra" 
of matrices in <p n to the length of the discretized non-orientable loop (f> n . We notice that, 
in each time step, the interaction process decreases the discretized loop length by the 
unit "2". The process which comes from the original action of matrix models extends 
the length of discretized loops. 

The definition of the F-P hamiltonian operator gives the precise definition of a field 
theory of second quantized non-orientable strings. In terms of the expectation value of 
an observable 0(0), a function of n 's, the F-P hamiltonian operator Hfp is defined by, 

< 0(O)|e-^O(0)|O >=< 0(^(r)) > 5 . (8) 

In R.H.S., 0|(t) denotes the solution of the Langevin equation (4) with the initial con- 
figuration 0(0) 7^ 0. The time evolution of R.H.S. is given by, 

< AO(0(r)) > f = < £ d m O{<Kr)) A0 m + W d„AO(0(r)) A0 m A0n > f +0(Ar 3 / 2 ) , 

m " m,n 

= -At < H fp (t)0(<P(t)) >t , 

(9) 
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where d n = By substituting the Langevin equation (4) and the noise correlation (5) 
into (9), we obtain 

H fp (t) = -J2 x nmr n , 

n>0 

1 1 ™~ 2 1 1 

N m 2 ^ 2 iV ^ 

(10) 

where n n = . To define the operator formalism corresponding to eq.(8), we introduce 
(f) m and 7r m as the creation and the annihilation operators for the loop with the length n, 
respectively. Then we assume the commutation relation [rt m , <f> n ] = 5 mn , and the existence 
of the vacuum, |0 >, with n m \0 >=< O|0 m = for m > . In the representation, 
< Q\ =< 0|e^m Qm7i " m and \Q >= U m 5((f) m - Q m )\0 > , the F-P hamiltonian operator 
Hfp in (8) is given by replacing <p rn — > <fi m , and n m — > n m in H FP in (10) with the same 
operator ordering. 

From the equality (8), the probability distribution function P((f), r) , which is defined 
by < O(0(r)) >^= /n n rf0 n O(0)P(0,r) , is given by, 

P(0,r)=< 0(O)|e-^|0> . (11) 

The initial distribution, P(0, 0) = U m 5(<f) m — m (O)) , represents the initial value of the 
solution of the Langevin equation (4). Eq.(ll) follows the Fokker-Planck equation for 
the probability distribution, 

AP(0, r) = +ArH FP P(<f ) , r) , (12) 
where H FP is the adjoint of H FP in (10), 
Hfp = n7r nX n , 

n>0 

1 1 n ~ 2 11 

Xn = ~ A?2 51 m7r m0m+n-2 + ~ X! 0fc0n-fc-2 + ~(™ ~ l)^T70n-2 
iv m Z fc=Q Z iV 

+ XI #cA+a • 

(13) 
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The remarkable observation is that it includes the Virasoro constraint Since the 
stochastic time evolution is generated by the noise essentially, the emergence of Virasoro 
constraint is traced to the noise correlations in eq. (5) which are equivalent to the insertion 
of matrices into the loop variable, M — > M + ArM m_1 , in <ft n . It generates the transfor- 
mation [— ArL m _ 2 ,0 n ] = nAr(j) m+n _2 , which corresponds to the noise correlation (5). 
In fact, for real symmetric matrix models ( non-orientable strings ), L n = —N 2 X n+ 2 
satisfies the Virasoro algebra without central extension, 

[L m , L n ] = (m - n)L m+n . (14) 

It is also worthwhile to note that the F-P equation (12) realizes the Virasoro constraint 
for the probability distribution. Namely, L n = N 2 X n+ 2 also satisfies the Virasoro algebra 
without central extension (14). Therefore, the F-P equation deduces a constraint equation 
for the distribution function even at the discretized level, justifying the generation of the 
partition function which satisfies the Virasoro constraint at the infinite stochastic time. 

L n fim P{(f>,r) = , for n = -1,0,1,... . (15) 

For hermitian matrix models, the Virasoro constraint for the partition function (15) was 
found as the S-D equation |i~E| . In the continuum limit, it deduces the continuum version 
of the Virasoro constraints jLTj]. The expressions (8) and (11) also give a constraint on 
possible initial condition dependence of the expectation value and the partition function 
at the infinite stochastic time limit, such as, lim^oo Hpp[(f)(0), gjhfi]Py>, t] = . This 
implies that these quantities may have the initial value dependence up to the solution of 
the Virasoro constraint. 



Continuum Limit and Continuum Field Theory of Non-Orientable Strings 

Now we take the continuum limit. First we introduce a length scale "a " and define 
the physical length of the loop created by <p n with I = na. Then we may redefine field 
variables and the stochastic time at the continuum limit as follows. 

G st = N- 2 a- D , 
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dr = a- 2+D / 2 Ar, 

m ee a-^, 

n(0 = a^V, 

(16) 

where we would like to keep the string coupling constant, G st , finite at the double scaling 



limit [IS]. For the existence of the smooth limit from the discretized stochastic time 



evolution to the "continuum" one, we require the condition, =j — 2 > 0. The scaling 
dimensions of all the quantities in (16) have been determined except the scaling dimension 
of the string coupling constant, D, by assuming [00, 

ArHpp = drTi.FP , 

(17) 

Then we obtain the continuum F-P hamiltonian, TLfp, from Hpp at the continuum limit, 

1 roo roo rl 

H non-ar. = / rf Z{2G st / dl'${l + l')l'U(l')m(l) + / dl'${l ~ l')$(l')lIL(l) 

2 Jo jo jo 

+ v /G^/$(/)/n(/)+p(/)n(/)} , 

(18) 

for non-orientable strings. By the redefinition (16), the F-P hamiltonian (18) is uniquely 
fixed at the continuum limit except the cosmological term. To specify the explicit form 
of the cosmological term p(l) in (18), we have to carefully evaluate the contribution 
which comes from the 3-point splitting interaction term and the matrix model potential, 
a l ~ D J2 a =o 9a4>n+a > at the double scaling limit of the real symmetric matrix model. Here 
we remember that the S-D equation for non-orientable string at large N limit takes 
exactly the same form as the orientable one under the suitable choice of the matrix 
model coupling constants. Since the continuum limit is taken by using the universal part 
of the disc amplitude 0, we naively expect p(l) takes the same form as that for orientable 
strings. 



To show explicitly this is indeed the case, we consider the simplest matrix potential 
given by, go = —1/2, g\ = g/2,g 2 = g% = ... = in (1), which corresponds to c = 0. Let 
us introduce the string field variable 

<j>{z) = £^>n = -U> ' 



N z- MN~V 2 



To take the continuum limit, we redefine the field variable, 

4>{z) = ^(z-gz 2 ) + CoZ^a 3/2 ^(u) 



(19) 



(20) 



where we have introduced the "renormalized" parameters, z = z c e au , and g = g c e~ Cia 1 , 

and the "continuum" stochastic time dr = co-2~ 2 a 1//2 Ar , where z c = 3^(3^ + 1) and 

q c = — ~~t" are the critical values and t denotes the cosmological constant. The constants 
y 2-3! 

Co and c\ are chosen for convenience. The scaling dimension of all the quantities have 
been determined so that the string coupling, 1/G st = clN 2 a 5 , is fixed at the double 



scaling limit [IS 



By using the Laplace transformation, 

POO 

= / dle- ul $(l) , (21) 
J 

in the Langevin equation derived for the redefined field variable $(w) in (20) with the 
same procedure as shown in (4), we obtain the following Langevin equation which is 
equivalent to the continuum F-P hamiltonian (18) |J, 

= d((l) + -dr{l [ l dl'm')${l - V) + p(l) + JG~ t l 2 <&(l)} , 
2 Jo v 

< dC(l)dC(l') > = 2drGJl' < $(Z + /')>> 



3/ 

P (i) = 3 a*(0- T *(i) 



(22) 
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for non-orientable string. It is consistent with the naive continuum limit of its discretized 
version (4) except the term p(l). As we have shown explicitly, the cosmological term p(l) 
takes the same form both for orientable and non-orientable strings. The field theory of 
non-orientable strings is also consistent with ref. 0] in transfer matrix formalism. The 
double scaling limit of the real symmetric matrix model has been studied in a quartic 
potential |13|, while our result shows that it happens in the cubic potential as well. We 



notice that the continuum F-P hamiltonian includes the continuum Virasoro generator 
£(/), 

roc I A 

c(i) = -{ / di , m+i')i i u{i , ) + —- di'$(i-i')<s>(i') 

JO Zbof JO 



. W(J) + -Lf) 



2V G s t 2G s t I 

(23) 

These generators satisfy the continuum Virasoro algebra, [£(/), C(l')\ — (I — l')C{l + V) . 
In the stochastic quantization view point, since the stochastic time scaling dimension is 
given by f -2 = | > for c = 0, we expect that the discretized version of the loop space 
Langevin equation for real symmetric matrix models may provide a possible method for 
numerical calculation of non-orientable 2D random surfaces to sum up the topologies of 
surfaces. In the next section, we extend the idea which we have learned in 2D case to 4D 
Euclidean quantum gravity. 



4D Quantum Gravity from Stochastic 3-Geometries 

Here we point out that the time evolution in 4D quantum gravity is described by a 

Langevin equation for 3-geometries in terms of the Ashtekar's canonical field variables 

by showing that the corresponding Fokker-Planck hamiltonian operator exactly recovers 

the hamiltonian of 4D Euclidean quantum gravity in the gauge N = 1 and N l = 00. 

The Hartle-Hawking type boundary condition is naturally imposed in this scheme by 

specifying the initial probability distribution functional. We also present the lattice reg- 
2 In 2D quantum gravity, it is pointed out that this gauge fixing recovers the time evolution in 
non-critical string field theories [?]. 
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ularization of this approach which defines a lattice regularization of Ashetekar's canonical 



formalism [14 



At first, We propose the basic Langevin equation for 3-geometry in terms of the 
Ashtekar's canonical variables |HJ to recover the hamiltonian of 4D quantum gravity with 
the corresponding F-P hamiltonian defined latter. The simplest form of the Langevin 
equation is defined by, 



A^(x,r) = AQ(x,r) 

K 

i 

2 



<AQ(x,T)A(»(y,T)> c = -Are abc <F^T)> ( S 3 (x-y) , (24) 



where A%(x, r) is a SU(2) gauge field in Euclidean Ashtekar's canonical formalism^. In 
this note, Latin indices "i,j,k,..." denote the spatial part of the spacetime coordinate in- 
dices. While the Latin letters "a,b,c,..." denote the spatial part of the internal indices. 
x, y, ... denote spatial spacetime coordinates. The one step time evolution is defined by 
AA$(x,t) = Af(x,T + At) — A*(x,t) in (24). The coupling constant k is defined by 
k = 16ttG with G, the gravitational ( Newton's ) constant in the natural unit % = c = 1. 
The noise variable in (24) is not a simple white noise. The expectation value of the 
R.H.S. of the noise correlation is understood to be taken with respect to the noises up to 
the one unit time step before, r — At, in the sense of Ito's stochastic calculus. It is also 
equivalent to require < AQ(t) >= in Ito's calculus. 

The invariant property of (24) is not apparent even if we introduce algebra-valued 1- 
form, A(x,t) = AfT a dx l ,A((x,t) = AQT a dx % , and algebra-valued curvature 2-form, 
F{x, t) = dA + AAA . Then (24) is rewritten by, 

AA(x,t) = A((x,t) , 
<A((x,r)AA((x,r)> = kAt5 3 (0)F(x, r) . (25) 

The R.H.S. of the noise correlation should be regularized in a gauge invariant way. The 
basic Langevin equation is manifestly covariant under the SU(2) local gauge transforma- 
tion, A^ix) — > Af(x) + Df b u b (x) . While it is not covariant under the spatial general 
3 We use the notation, F£ = - d.Af + e abc A b 4 A° 4 , and Df c = 5 ac di + e abc A b . The field variable 
Af(x) is real in the Euclidean Ashtekar's formalism. 
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coordinate transformation, Af(x) — > Af(x) + F^{x) , due to the appearance of the 
divergent term, kAtS 3 (0)F^ : ' , in the transformation of the R.H.S of (24). It is formally 
cancelled by adding a term in the R.H.S. of the Langevin equation (24) which comes from 
the invariant path-integral measure we discuss latter. 

In terms of the solution of the Langevin equation, the following equality holds. 

<U x>i> J(At c (x,r) - Af inal \x))> c 

^ j^final)\ e -TH FP [Tt,A}\ ^{initial) ^ 

= < j\( initial ) \q- tH fp[AM tj^ifinal) ^ ^g) 

In the L.H.S., A^(x, r) denotes the solution of the Langevin equation with the initial 
condition, A?(x,0) = Af nitial \x) . In the R.H.S., H fp [tt,A] and H FP [A, ft] are defined 
by, 

H F p[*,A) = JSxe abc K^)^)F^x) , 

H FP [A,fr] = -^Jd 3 xe abc F^x)r a (x)7Tl(x). (27) 

To show the equality (26), the commutation relation, 

K(x) ,A b (y)]=S b a S]S%x-y) , (28) 

and the vacuum, |0 > with ft l a \0 >=< 0\A b = 0, have been assumed ( see also eq.(8)) |J. 
The Fokker-Planck hamiltonians (27) are just the hamiltonian for 4D quantum gravity 



without cosmological term in Ashtekar's variables |19] with different operator orderings. 
The operator 7r* is interpreted as the "triad "in the Euclidean Ashtekar's formalism. 
NamelyQ, fc J a (x) = ^e^(x) . Stochastic time evolution with the Langevin equation (24) 
corresponds to the Euclidean time evolution in 4D quantum gravity with gauge fixing, 
N — 1 and N l = 0. In the stochastic process (24), the "triad "plays the role of the noise 
variable and the scale of the fluctuation is characterized by the curvature. Though the two 
expressions in (26) are precisely equivalent, the second definition of the F-P hamiltonian 



4 e l a = q 1 / 2 e l a , where q = det((to) with spatial metric q l i — e l a e 3 a 
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in (27) has an advantage for further discussion. This is because the vacuum satisfies the 
( local ) hamiltonian constraint with the operator ordering in Hpp, < Q\H(A(x), 7t(x)) = 
H(A(x),tt(x))\0 >= , where H FP [A,n] = Jd 3 xH(A(x), n(x)) . It should be noticed 
that the operator ordering of the F-P hamiltonian operator H FP is deferent from that 
appeared in the F-P equation. 

Let us consider the initial distribution dependence of the probability distribution 
functional by averaging the expectation value (26) with respect to the initial probabil- 
ity distribution. It is defined by integrating out the initial configuration A^ mitml ' > (x) , 
on which the solution of the Langevin equation A^(x,t) depends, with the distribu- 
tion, p[A( imt%al \0], in the L.H.S. of (26). It gives a generalized form of the distribution 
functional P[A,r], which is defined by < 0[A c ](t) > c = JVAO[A]P[A,t], as follows 

P[A, t]= J VA imtial < A initial \P[A, ti\e~ THFp[A ^\A > . (29) 

For an arbitrary observable 0[A], the average with respect to the initial value distribution 
also gives, 



< 



0[A c (t)] > c = JvA initial < A initial \P[A,0]e- THFp[A '* ] O[A}\0 > . (30) 



In the definition of the expectation value in the L.H.S. of (30), the average is also taken 
with respect to the initial values with the distribution p[A mitml , 0]. For example, eq.(26) 
is given for 0[A] = 11^5 (X a (x) - Af final \x)), with the initial distribution, P[A,0] = 
Rx,a,i$(Af (x) — A^ miUal \xfj . From (30), the time evolution equation for the expectation 
values of observables is given by, 

j- < 0[A c (x,r)] > c =< H FP [A%x), j^PiA^x)]^^ > c . (31) 

The initial condition dependence in the amplitude (30) leads a restriction on the 
boundary condition with respect to the Euclidean time. We notice that eq.(30) also gives 
a constraint for initial distribution, 



^- < 0[A c {x,t)] > c = [vA in P[A in , 0}H FP [A in , — -} < A in \e- THFp[A ^ ] 0[A]\A > .(32) 

(XT J 1\ 
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The existence of the equilibrium limit requires that the R.H.S. in (32) should be zero at 
the infinite stochastic time. A trivial solution of this constraint is the vanishing curvature 
at any spatial points, P[A, 0] = H x ,a,ijfi \F^(x)J . This initial distribution, however, should 
be excluded because the R.H.S. of (31) is identically zero even at finite stochastic time 
and there is no time development. In general, if we choose a solution of the hamiltonian 
constraint as the initial value distribution, obviously there is no time evolution in the 
corresponding F-P equation. Especially, the constraint does not allow us to solve the 
Langevin equation with the initial condition mfaa ^(a;) = 0. 

To specify the physical boundary condition and define a class of solutions for the 
local hamiltonian constraint of 4D quantum gravity in this context, we may choose the 
following initial condition which generates a nontrivial time evolution as an analogue in 
2D case, 



The spatial coordinates, x = zq, is identified to the point where the 3D manifold is 
absorbed into nothing in the sense of Hartle-Hawking type boundary condition. The local 
hamiltonian constraint is broken at this point but it may be recovered at an equilibrium 
limit. 

Let us consider the gauge invariant path-integral measure in the Ashtekar's variable, 
which are implicitly assumed in the expectation value (30). One way to specify the 
path- integral measure is to introduce a regularization for the noise correlation (24). In 
the following, we show that an extra term is necessary in a lattice regularization of the 
Langevin equation (24) to define an invariant and well-defined measure by identifying 
the noise correlation in (24) to a "superspace "metric with Ashtekar's variable. 

By using the lattice regularization of the Langevin equation and the noise correlation 
(24) , the invariant property of the Langevin equation in the sense of Ito's calculus 
naturally introduces the path-integral measure in the configuration space with Ashtekar 
variables. It is given by $M 



AU(x,t) al3 (r) = AW(x,tU(r) + ArlGl 1 / 2 (\G\~ 1/2 G[x, i\ V, jW) • (34) 




(33) 
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The regularized Langevin equation describes the one step time evolution of the link 
variables^, U(x, z)(r + At) = U(x,i)(t) + AU(x,i)(t) . The dynamical variable and 
the noise variable, U(x,i) and AW(x, i) respectively, have been assigned on the link of 
3-dimensional lattice, which is specified by the site x and its nearest neighbor in the 
i-th direction denoted by x + i. U(x, i) is an element of SU(2) group in the adjoint 
representation, while noise AW(x, i) is algebra valued. The quantity Q[x, i; y,j] a p- tl s , is 
interpreted as the inverse of the "superspace "metric and \Q\ denotes its determinant. The 
superspace is spanned by the configuration {U a/ 3(x, i)}. The inverse of the superspace 
metric is given in the following regularized noise correlation. 

< AW(x, i) al3 (r)AW(y, j) l8 (r) >= At < g[x, i; y, ]] a p- n s{r) > , (35) 

where, 

G[x,i;yJ}af3;~fS = S x+ i y+ 'j5 0S {U(x,i)U(x + i,-j) -U(x,-j)U(x-j,i)} cn 

+ 5 x+i,y 6 Pl{ U ( X ^) U ( X -U(x,j)U(x+J,i)} a6 

+ 5 Xi y + j5 a s{U (x + i, -i)U(x, -j) -U(x + i, -j)U(x + -i)}^ 
+ 6 x> y6 ai {U(x + i,-i)U(x,j)-U(x + i,j)U(x + i+j,-i)} f3S . (36) 

The regularized Langevin equation (34) is invariant under the "general coordinate trans- 
formation "in superspace, {U(x,i)}. 

U{x,i) a/3 ^V{x,i) af3 [U] , (37) 

where V a p(x,i) is also an element of SU(2) group in the adjoint representation and an 
arbitrary functional of U^siy, j)- The second term in the R.H.S. in (34) is necessary for the 
invariance in Ito's calculus |14j]. The role of the contribution from the invariant measure 
was first clarified in the case of the Langevin equation for a particle moving on a Riemann 
surface |23|] ( see also Ref. J25| and the Appendix in Ref. |26| ). The similar argument 



5 For the detailed discussion to keep the dynamical variable within the element of SU(2) group along 
the time development and for the derivation of the corresponding lattice regularized F-P hamiltonian, 
see Ref. f3 
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is formally possible even in the non-regularized version (24) if we identify the R.H.S. 
of noise correlation in (24) to the continuum superspace metric. Therefore, the naive 
continuum limit of these equations, (34) (35) and (36), coincides with those in eq.(24) 
except a divergent term with 5 3 (0) which comes from the second term in the R.H.S. of the 
Langevin equation (34). Though the divergent term actually represents the contribution 
from the path-integral measure and it is necessary for the invariant property of the 
Langevin equation in a formal sense, it is not well-defined without regularization. Thus 
the regularized Langevin equation (34) provides a possible basis for numerical simulation 
on 4D quantum gravity. 

We comment on the two important consequences of the lattice regularized Langevin 
equation |T4[. One is the corresponding equilibrium distribution, 



lim P[U,r] = g~ 1 ' 2 , (38) 

r— >oo 

which defines a measure invariant under the superspace general coordinate transformation 
(37), 

VU = g-y 2 U Unk dU , (39) 

in the regularized version of the expectation values (29) and (30). The other is the 
Schwinger-Dyson equation in this context. It is given by, 

< \G\ l/% WTT ~\ (|gr 1/2 £M^W) >=0 , (40) 



at an equilibrium limit. 



Conclusion 

In conclusion, we have shown that the Langevin equation for real symmetric matrix 
models written by the loop variables defines the time evolution of non-orientable strings 
which defines non-orientable 2D surfaces at both discretized and continuum levels. The 
partition function in loop space satisfies the Virasoro constraint at the equilibrium limit 
in both discretized and continuum level. 
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The idea has been extended to 4D Euclidean quantum gravity. In the Langevin 
equation (24) and the regularized one (34), the "triad "in Ashtekar's variables are re- 
alized as noise variables, which presumably represents the stochastic 3-geometries. The 
corresponding F-P hamiltonian is equivalent to the hamiltonian of 4D Euclidean quan- 
tum gravity in the temporal gauge. It allows us to interpret the stochastic time as the 
Euclidean proper time. The strategy we would like to adopt here is to characterize the 3- 
boundaries in the 4D spacetime by using the solution for both momentum constraint and 
Gauss law constraint as the initial distribution and the observables. As it is clear from 
(30), initial distributions and observables are the key quantities to specify 3-boundaries on 
a 4D spacetime in this approach. The method developed in 2D quantum gravity in loop 
space indicates that the expectation values of such quantities satisfy all the constraints 
in 4D quantum gravity at the equilibrium limit. There are some candidates useful to 
characterize 3-boundaries, such as the extrinsic curvature term, 3D Chern-Simons term, 
topological invariants and loop variables for which the present formalism can 
be applied. For these observables, the gauge fixing terms for the spatial general co- 
ordinate invariance and local Lorentz invariance, which are introduced as a drift force 
in the Langevin equation (24) following a standard method []2"3| ||25[| , do not change the 



expectation value of observables at the equilibrium limit fL4 |. 

Though there is no well-defined drift force in the Langevin equation (24) without 
regularization of path-integral measure, there is another machinery to introduce a drift 
force. In the Langevin equation for observables which is equivalent to the F-P equation 
for the expectation value of observables (30), a drift force appears effectively as a direct 
consequence of the fact that Ito's stochastic calculous pick up the Jacobian factor which 
comes from the change of variables from the gauge fields to observables^ ||][]25|. In the 
6 The simplest example is 2D Euclidean Yang- Mills theory. It can be defined by the same Langevin 
equation for SU(N) as (24) with 1-dimensional field variables A a (x,T) and the white noise correlation, 
< A( a (x, r)A( b (y, r) >£= J? AT5 ab 5(x— y) . The noise is translated as the canonical momentum variable. 
Then the Langevin equation for the expectation value of Wilson loops defines the collective field theory 
for 2D Yang-Mills field. 
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lattice regularized Langevin equation (34), the contribution from the invariant measure 
introduces a well-defined drift force. 

The problem of the present scheme is that the noise correlation in the basic Langevin 
equation (24) is not positive definite. It would force the Langevin equation (24) to be 
complex, although the field variables are real in the Euclidean Ashtekar's formalism. The 
point would be main difficulty for the numerical analysis in this scheme. One way to deal 
with the problem may be to extend the Langevin equation (24) to a class of more general 
gauge fixing. It is always possible by multiplying the noise correlation in (24) with lapse 
function. Then one may chose the lapse function so that the noise correlation keeps the 
value to be positive definite. It is an open question if this gauge fixing procedure, a choice 
of non-trivial lapse function, make sense in numerical simulation. Another way may be 
to consider the Wick rotation of the conformal mode. Apart from these questions, the 
description with the Langevin equation has a topological feature in the sense of Nicoli- 
Langevin map. Such a topological feature would relax another difficulty, renormalizability 
of quantum gravity. We hope that the approach is useful for deeper understanding of 
quantum gravity. 
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